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Abstract 

The question of whether hydrogen atoms can exist or not in spaces with 
a number of dimensions (D) greater than 3 is revisited. The lowest quan- 
tum mechanical stable states and the corresponding wave functions are de- 
termined by applying Numerov's method to solve Schrodinger's equation. 
States for different angular momentum quantum number and dimensional- 
ity are considered. One is lead to the result that hydrogen atoms in higher 
dimensions could actually exist. The most probable distance between the 
electron and the nucleus are then computed as a function of D showing the 
possibility of tiny confined states. 
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1. Introduction 

Unification efforts based on superstring theories reinforced the idea that 
physical space dimensionality should be somehow explained in the framework 
of theoretical physicSjCven if some dimensions are supposed to curled up into 
a compact manifold [li. A modern and comprehensive survey of dimension- 
ality can be found in [2] . Searches for experimental evidences of extra dimen- 
sions are part of contemporary tendencies for investigating physics "beyond 
the Standard Model" in Collider Physics jsf. 

The idea that a particular physical law should depend on space dimen- 
sionality can be traced back to a philosophical speculation due to the young 
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Kant namely, the suggestion that the Newtonian gravitational force de- 
pends on the three dimensionality of space. This insight was indeed prolif- 
erous and had inspired, for example, Ehrenfest, who was the first to give a 
mathematical ground to Kant's hint by discussing the mechanical stability 
of the formal solutions of the planetary motion in D dimensions, assuming 
it is still described by a Poisson equation generalized to , as well as the 
implications of higher dimensions on Bohr's atomic model js], [oj. 

Those papers bring the belief that one can learn about space dimen- 
sionality from a class of generalized physical equations. In particular, this 
supposition was the starting point of an article by Tangherlini [7], where 
the problem of Schrodinger's hydrogen atom in D dimensions was formally 
treated for the first time. Although the basic assumption that one can prove 
the three-fold nature of space from the generalized Schrodinger's equation 
was shown to have epistemological limitations jsj, the general ideas that it 
is possible to understand how the structure of a particular physical equation 
depends on dimensionality, or how dimensionality itself can depend on the 
spatial scale or even how a physical phenomenon could vary by changing 
the number of dimensions or other topological feature of space are still of 
considerable interest. 

In a nutshell, the study of such a relationship between the structure of 
physical laws and space dimensionality is particularly attractive in the frame- 
work of the general theoretical scenario where extra space-time dimensions 
should play an important role on several attempts to unify the fundamen- 
tal forces. Also the expectation from various scenarios of quantum gravity 
that the space-time dimension seems to rely on the size of the probed region, 
being somewhat smaller than four at small scales and monotonically raises 
with increasing the size of region ^ can be argued to justify the investigation 
on how physical laws depend on or are entangled in space-time dimension- 
ality. As a few examples of such a kind of dependence, one can quote the 
Casimir Effect [ill . [12I , the existence of scalar field lumps [lij and how 



recent data on cosmic microwave background can be used to settle an upper 
limit for fractal space dimensions [l^. In any case, it could be possible to 
learn more about physics in 3 dimensions by investigating its generalization 
to spaces of higher dimensions. 
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2. Former predictions 



Are there stable hydrogen atoms in flat spaces if the number of dimensions 
is greater than three? Both positive and negative answers are found in the 
literature. One can try to summarize the main results as follows. 

In his semi-classical approach, Ehrenfest, based on Bohr quantization of 
circular atomic orbits, characterized by the total quantum number n, had 
shown that for orbits in D dimensions, the Bohr result for D = 3, 

1 m 

should be replaced (for D > 2) hj 



(ri\ {D — A) 2D-4 (ra\%^ , . 

Elf^ = ^ — -Lxn-^x( — ) xe^ (1) 



2{D-2) 

So, according to this result, the energy becomes positive for D > 5 and 
increases to infinity while the orbits draw closer and closer to the nucleus. 

To the best of our knowledge, Tangherlini was the first, in 1963, to for- 
mally treat the problem of the hydrogen atom from the point of view of wave 
mechanics 0. The problem was considered by him from the point of view 
of stability arguments and he did not solve the Schrodinger equation for an 
electrostatic potential energy of the type V = e/{D - 2)r^-2. The author 
limited himself to give a semi-classical argument based on the mechanical 
stability of the classical orbits. He immediately claimed that it is clear that 
for the cases D > 5, the energy levels have a point of accumulation at minus 
infinity: i.e., r = is not as regular point; and hence there are no stable 
bound states. The case D = 4, can also be excluded by standard arguments. 



The same argument is sustained in Ref. [15|]. Nevertheless, to the best of 
our knowledge, this affirmative answer should be taken with reserve. Indeed, 
from the general theory of differential equations, it is well known that when 
in the equation 

y" + P{x)y' + Q{x)y = 

one or both coefficient functions P{x) and Q{x) fails to be analytic at Xo, this 
point is a singular point, and it is said to be irregular so far {x — Xo)P{x) 
and (x — Xo^Qix) are not analytic. Whenever this happens, one can safely 
conclude that there might not exist a Frobenius series solution for the dif- 
ferential equation [l6j, and that is all. In fact, a counterexample can be 
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found in 17|, |l8| where analytic study of the dominant effects of a polariza- 
tion potential that depends on r (distance between an incident charge and 
a multi-electron atom) like 1 / r'^ was done. Corrections to the inverse sixth- 



power term were treated in [19[. Indeed, it was shown that the Schrodinger 
equation for inverse fourth- and sixth-power potentials reduces to peculiar 
cases of the double-confluent Heun equation and its Ince's limit, respectively 
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In 1971, the same problem was investigated supporting the conclusion 
against the existence of stable atoms and atomic structures for D > 3 with 



negative energy spectrum, although positive spectrum was not analyzed [21 



In 2005, it was claimed that for one electron atoms in dimensions greater than 



three there is no normalized wave function corresponding to bound states [22 

However, Burgbacher and collaborators, in 1999, have shown that it is 
possible to have stable hydrogen atoms in higher dimensional non-compact 
They admitted that the kinetic energy in the Schrodinger oper- 
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spaces 

ator has the usual form described by the D -dimensional Laplacian and that 
the electrostatic interaction in the Schrodinger equation has the same form 
irrespective of the spatial dimensions. Instead, they are obliged to change 
the structure of Maxwell equations in higher dimensional spaces. In anyway, 
so far Maxwell's theory is concerned, one should not forget Weyl's classical 
result which showed that only in a (3+1) dimensional space-times can such 
electromagnetic theory be derived from a simple gauge invariant action inte- 



gral, having a Lagrangean density which is conformally invariant [2^, |25|, |26 



On the other hand, exact solutions for the radial Schrodinger equation in 



D dirnensions were found for potential of the type V{r) 



ar 



'^+br '^+cr ^ + 



dr-^ |271]. In the paper [2^ 



the discrete spectrum of the Schrodinger equation 
corresponding to the potential V{r) ~ 1/r, is treated in D dimensions. The 
authors considered also in details the Cornell potential, V{r) = —a/r + br. 
Algebraic solutions for the supersymmetric hydrogen atom with the potential 
V{r) ~ 1/r are given in 29 . 

Another sort of approach considers the possibility of confined hydrogen 



atoms in higher space dimensions by spherical cavities [30|, |3l| • The ground 
state energies were computed in impenetrable spherical cavities for different 
dimensions by using the Coulombian potential ~ 1/r |31]. 

All those ideas have motivated several approaches to the problem but 
with conflicting conclusions. In many of them, the kinetic part of Schrodinger 
operator is generalized while the Coulombian energy potential is taken as in 
the case of three dimensions. 
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On the other side, in Ref. [32[, it was shown that the Dirac equation in D 
dimensions, generahzed for a Coulombian potential (which is a superposition 
of a Lorentz scalar term, V{r) = —a/r, and the temporal component of 
a Lorentz quadri- vector of the type V{r) = —b/r), has an exact solution. 
There are also a number of works that generalize Dirac equation for one 
electron atoms under the action of an attractive central potential. In many 
of them, the Coulombian potential form in three dimensions, V{r) ^ 1/r, is 



maintained no matter the value of D 33 , 34 , 35 . Those works can also be 



regarded as a further motivation to revisit the problem of the existence (or 
not) of hydrogen atoms in higher dimension spaces. Since in this Letter we are 
not considering external magnetic fields and since spin effects on the energy 
spectrum are known to be very small (of the order of ~ 10~^) in the Pauli 
equation, this spectrum can be determined in a very good approximation by 
Schrodinger equation. 

Admitting the hydrogen atom to be described by a generalized Schrodinger 
equation, is it true that whenever one postulates the validity of Gauss law 
in D dimension non-compact spaces - which requires the replacement of the 
potential energy by a function that goes like ~ l/r^~^, in addition to the 
generalization of the kinetic part of Schrodinger's Hamiltonian operator -, 
the very existence of hydrogen atoms or their stability is lost? This is what 
will be investigated in the next Section. 



3. General results 

In a flat (Euclidean) space, the radial Schrodinger equation in D dimen- 
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sions can be inferred from the D-dimensional Laplacian operator V?r)-, in 



spherical coordinates, which is given by 

where r is the radial coordinate, and Vt^ is the Laplacian operator on the unit 
hyper-sphere R^~^. The solution of the generalized Schrodinger equation, 
ip{f)., can be expressed in terms oi D — 1 angles as 



^(r) = i?(r)F,(e,0i,02,---0 



D-2 J 



The functions are the Gegenbauer polynomials, which are the angular 
solutions of the equation 
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while the radial solution R{r) satisfies 



d'^Rh) D-ldR 

— H h 

dr^ r dr 

,2 



2m 



i{£ + D- 2) 



2m r 



2 



i? = (2) 



where C/(d) = -e(D) V(£,)(r). 

Making the choice analogous to the three dimensional case, 

u{r) = r(^-^)/'i?(r) (3) 

the first order derivative can be eliminated, giving rise to 

^ + |^P-Mr))]«(r) = (4) 

where the effective potential energy is given by the electrostatic energy plus 
the generalized centrifugal term 



with m and e(D) being, respectively, the mass and the charge of the elec- 
tron and the values of j depend on the orbital angular momentum quantum 
number / and on the number of dimensions D as 

+ ^ (6) 

Up to this point there is a consensus, but the choice of V{r) is not unique 
in the literature. Often, the usual Coulombian potential established for 
D — 3 is assumed to be valid for an arbitrary D. Instead, we will use 
throughout this Letter the formula that guarantees the electric charge con- 
servation in a D dimensional space, or, equivalently, the generalized Poisson 
equation, where the power of the Laplacian operator is maintained and just 
the number of independent spatial coordinates varies, i.e., 



2T{D/2) ep) 
7r(-D-2)/2 (D_2)r^-2 



Once the potential is fixed, it is straightforward to see that it does not 
belong to the class of potential to which the Schrodinger equation is exactly 
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solved 361]. Therefore, at this point, Eq. will be solved for the potential 



energy given by and ([7]) by applying the Numerov numerical method 



371 . 138l . |39[ to solve an eigenvalue differential equation. A specific program 
was developed by the authors in C + + language and both calculations and 
graphics were done by using the CERN/ROOT package. 

In this Letter, the Schrodinger equation is solved numerically for few 
values of angular momentum, i.e., £ = 0,1, 2, 3, allowing space dimensionality 
to vary between the interval 5 < D < 10. 

Before presenting the numerical results, one might say that, so far con- 
cerning the problem of space dimensionality, one should adopt something 
similar to the cosmological principle, according to which one should expect 
that the laws of physics, as determined in our neighborhood, are valid in all 
regions of the Universe and in all moments of its history, independently of 
the space-time scale which is being probed. An analogous hypothesis should 
be assumed here as in all literature, namely, that the physical law established 
in D = 3 will be valid for different values of the dimensionality D and that 
the numerical values of the physical constants should not vary significantly. 
Otherwise it will be impossible to get any numerical prediction. 

The present calculations depend on the numerical value of the electric 
charge, e(£)), in D dimensions. It will be assumed throughout this Letter 
that e(£i) has the same value of e measured in three dimensions. This should 
be justified from the results of ref. jioj, where it is shown that the numerical 
value of the generalized fine structure function for D dimensions is very close 
to the three dimensional one, 1/137, and also Planck constant h^o) did not 
vary significantly with space dimensionality. 

With this assumption, the effective potential energy is given by fig. [Hand 
the energy eigenvalues are calculated for the first four principal quantum 
numbers {n = 1,2,3,4) and the lowest orbital angular momentum i = 0,1, 
varying D from 5 to 10. The results are shown in tabled] and table [2l 

How these lowest energy eigenvalues (corresponding to the states i = 1) 
grow with the increasing of the space dimensionality is shown in fig. [2l 

It is important to stress that all those values are very small compared 
to the maximum value of the effective potential (Um) for each value of D. 
Indeed, considering the i = 1 state, the ratio Ei/Um goes from 0.1% to 0.6% 
when D goes from 5 to 10. These results suggest a very low probability of 
tunnehng effect. 

In addition, we studied how the energy spectrum depends on i, for a 
particular D. For example, for D = Q and £ = 0,1,2,3, the energy eigenvalues 
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Table 1: Energy eigenvalues corresponding to principal quantum numbers n = 1,2,3,4 
and ^ = 0, for dimensions 5 < < 10. 



En [eV] 


D = 5 


6 


7 


8 


9 


10 


El 


0.31 


0.41 


0.51 


0.62 


0.74 


0.88 


E2 


1.77 


2.05 


2.30 


2.57 


2.86 


3.14 


E3 


4.39 


4.86 


5.30 


5.72 


6.12 


6.56 


Ei 


8.16 


8.88 


9.49 


10.03 


10.63 


11.20 



Table 2: Energy eigenvalues corresponding to principal quantum numbers n = 1,2,3,4 
and ^ = 1, for dimensions 5 < < 10. 



En [eV] 


D = 5 


6 


7 


8 


9 


10 


El 


0.51 


0.61 


0.76 


0.89 


1.04 


1.20 


E2 


2.28 


2.58 


2.85 


3.14 


3.42 


3.77 


E3 


5.25 


5.70 


6.13 


6.57 


7.03 


7.47 


Ei 


9.40 


10.06 


10.63 


11.21 


11.80 


12.51 



are shown in table El 

All the regular radial wave functions were determined and just few of 
them are plotted (with arbitrary normalizations) in Fig. |3] (£ = 0), Fig. H] 
(£ = 1) and Fig. |5] = 3), all of them for the particular value D = 6. 

Once the radial wave functions, the mean value of the radius of the elec- 
tron (r) for a particular dimensionality D, given by 

(r) = j rRj{r)r^-^ dr dQD-i (8) 

was numerically calculated by solving the above integral by using the rejec- 
tion Monte Carlo Method. 

Just the predicted values for (r) and the respective errors (for the state £ = 
1) as a function of D are shown in table HI where the values are expressed in 
terms of the numerical value of the Bohr radius qb, equal to 0.529 x 10~* cm. 
It should be noticed that for the quantity {r) /as is multiplied by a factor 5 
going from D = 5 to D = 10. 

In fig. [HI it is plotted the dependence of (r) (normalized to the Bohr radius 
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Table 3: Energy eigenvalues corresponding to principal quantum numbers n = 1,2,3,4 
and^ = 0,l,2,3, for £) = 6. 



En [eV] 


L = 


L = 1 


L = 2 


L = 3 


El 


0.410 


0.615 


0.976 


1.085 


E2 


2.051 


2.582 


3.183 


3.406 


E3 


4.864 


5.702 


6.664 


6.971 


E4 


8.878 


10.057 


11.235 


11.666 



Table 4: Dependence of (r) upon the dimension D, for ^ = 1 states. Errors for each 
prediction are shown in the last column. 



D 






5 


0.045 


0.001 


6 


0.090 


0.002 


7 


0.120 


0.002 


8 


0.174 


0.003 


9 


0.223 


0.003 


10 


0.252 


0.003 



m. D — ?)) on D in the case I —1. 

The mean value of the distance between the electron and the proton in 
the hydrogen atom tends to grow by increasing the number of dimensions for 
-D > 5 (for D = 4: there is no state). This means that the distance between 
the electron and the proton, for = 5, is an order of magnitude greater 
than in the case D — Notice that this result is qualitatively different to 
Ehrenfest's statement that the orbits draw closer and closer to the nucleus 
as the value of D increases. 

4. Final comments 

In this Letter it was shown that there is no state of a hydrogen atom 
with negative eigenvalues when space dimensionality is higher than 3. It is 
important to point out that for D = 4 there are no solutions for this kind 
of atom. Additionally, for dimensions D > 5, it is argued that there are 



9 



stable states with positive energy, described by well behaved wave functions, 
corresponding to an electron confined in a potential created by the central 
proton. The first energy levels are shown to be very small compared to the 
barrier high in order the electron could be confined in the proton potential, 
and the mean distance between the electron and the nucleus is about ten 
times greater than the correspondent distance in three dimensions. 

The amount of energy necessary to ionize the atom in higher dimensions 
may be numerically estimated by our method, by computing all possible 
stable states with energies less than the maximum value of the confined 
potential generated by the central proton. For instance, see tabled for D = 
5, the value of i?4 is indeed the biggest one still less than the maximum value 
of the potential. Thus, in this case, one can estimate the ionization energy as 
something close to this value. For other values of space dimensionality, more 
states fulfilling this condition can be found, even if they were not shown in 
table m 

In summary, one can ask what could be learned from the results pre- 
sented here. Admitting the hydrogen atom is described by the generalized 
Schrodinger equation no matter what space dimensionality is, it was found 
that now there are stable mathematical solutions not only for D = 3, against 
to what was claimed by different authors in the past. Indeed, it was shown 
that there are stable states even for higher dimensional Euclidean spaces, 
but all with positive energies. Therefore, the experimental observation that 
energy of the bound state hydrogen atom is of the order of —13.6 eV seems 
to indicate that nature should somehow prefer tridimensional space. 
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Figure 1: Potential energy of the £ = 1 state of a hydrogen atom as a function of r/as for 
the space dimensionahty 3 < £> < 10. 
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Dimension 



Figure 2: Energy eigenvalues for £ ~ 1 state of an hydrogen atom as a function of space 
dimensionality. The black squares refer to the Ei eigenvalues; the red triangles to the E2; 
the green, to E3 and the blue open circles to E4. 
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r/a^ 

Figure 3: Radial wave function i?f=o for the first four states of the hydrogen atom in the 
case D = 6. The color legend is: black curve refers to Ei eigenvalue; red to E2; blue to 
E3 and green to E4. 
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Figure 4: Radial wave function Ri=i for the first four states of the hydrogen atom in the 
case D — 6. The color legend is: black curve refers to Ei eigenvalue; red to E2; green to 
E3 and blue to E4. 
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Figure 5: Radial wave function Ri=3 for the first four states of the hydrogen atom in the 
case D — 6. 
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Space Dimension 



Figure 6: Dependence of the mean values of the ground state 
atom radius upon space dimensionality. 
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